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Abstract 

We study the band-centre anomaly in the one-dimensional Ander¬ 
son model with weak correlated disorder. Our analysis is based on the 
Hamiltonian map approach; the correspondence between the discrete 
model and its continuous counterpart is discussed in detail. We obtain 
analytical expressions of the localisation length and of the invariant 
measure of the phase variable, valid for energies in a neighbourhood of 
the band centre. By applying these general results to specific forms of 
correlated disorder, we show how correlations can enhance or suppress 
the anomaly at the band centre. 

Pacs: 71.23.An, 72.15.Rn, 05.40.-a 


1 Introduction 

The one-dimensional (ID) Anderson model, first introdnced in 1958 [1], re¬ 
mains the focus of active research. The analysis of this model showed that 
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all eigenstates are exponentially localised for infinitely large samples (see, 
e.g., the reviews BE]). For the case of weak and uncorrelated disorder, 
Thouless derived the expression for the localisation length that now bears 
his name |lj. Numerical calculations, however, revealed that Thouless’ ex¬ 
pression fails at the band centre, where the actual localisation length was 
numerically found to be different from the expected value [3]. As was un¬ 
derstood by Kappus and Wegner [B], this discrepancy is due to the failure 
of the Born approximation; with the use of a quite sophisticated method, 
Kappus and Wegner derived an approximate expression which explained the 
numerical data reported in Ref. [5]. Later on, using a different approach, 
Derrida and Gardner showed [7] that this anomaly is a resonance effect that 
emerges at the band centre, = 0. They also found another anomaly for 
£' = ±1 and suggested that similar anomalies should appear for other reson¬ 
ant energies dehned hy E = 2 cos{'n'p/q), where p and q are integer number. 
The results of Ref. |7] were subsequently extended to a whole neighbourhood 
of the band centre mi (see also the discussion in Ref. [TO]). 

The interest in the anomalies of the Anderson model was raised signi- 
hcantly after these anomalies were linked to the so-called single-parameter 
scaling theory (SPS). The SPS theory was proposed in Ref. [HI [12] and 
was based on the random phase approximation (RPA) for the phases of the 
scattered waves. In the ID case, the SPS theory can be reduced to the 
statement that all moments of the transmission coefficient through a random 
barrier can be expressed in terms of the first two moments. As a result, the 
whole distribution of the transmission coefficient is expected to depend on 
the ratio between the localisation length and the size of the random chain. 
However, it was numerically shown long ago |T3| that the RPA fails at the 
band centre of the Anderson model. Thus, the band centre anomaly (as well 
as the other resonances) can serve as a touchstone for the study of the SPS 
hypothesis. As was recently shown in Ref. [T3|, the SPS hypothesis is viol¬ 
ated also for the resonance at £' = ±1 (see also the discussion in Refs. na). 
Since the phase distribution is not flat whenever the energy lies in a whole 
neighbourhood of the resonant values, one can expect that the SPS hypo¬ 
thesis should not be rigorously valid in the whole energy band of the ID 
Anderson model with weak disorder. This “anti-SPS” hypothesis is still not 
proved; the study of the anomalies of the Anderson model for E = 0 and 
E = ±1 can provide a way to test the validity of this conjecture. 

The works on the band-centre anomaly mentioned so far considered the 
Anderson model with uncorrelated disorder. In the literature, almost no 
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attention has been given to the case of correlated disorder, the main exception 
being Ref. [IB] which, however, was focused on specihc cases of very short- 
ranged correlations. Only very recently the effects of correlated disorder on 
the localisation of the band-centre states have begun to be investigated, in the 
wake of the discovery that the band-centre anomaly can be strongly reduced 
if the site energies exhibit exponentially decaying, positive correlations na. 

In this work we derive analytical results that show how the band-centre 
anomaly is modihed when disorder possesses arbitrary correlations. The gen¬ 
eral expressions that we obtained allow us to treat the correlations considered 
in [a as a particular case and to show how, in general, different kinds of dis¬ 
order correlations can either enhance or suppress the anomaly at the centre 
of the energy band. To derive an analytical expression for the localisation 
length in a neighbourhood of the band centre, we relied on the Hamiltonian 
map formalism [181 [10] and we replaced a map for the angle variable with 
its continuum limit. In order to do so, we had to establish a rigorous cor¬ 
respondence between random maps with intercorrelated coloured noises and 
stochastic differential equations with the same features. This led us to derive 
a discrete integration scheme for stochastic differential equations with self- 
and cross-correlated noises. This scheme is another important result of the 
present paper. 

To analyse the band-centre anomaly we considered the behaviour of the 
localisation length and that of the probability distribution for the angle vari¬ 
able (which is equivalent to the scattering phase). Our analytical results show 
that disorder correlations shape the localisation length in a twofold way: on 
the one hand, the modulation of the phase distribution causes a deviation 
of the actual Lyapunov exponent from the value predicted by the expres¬ 
sion hrst derived by Izrailev and Krokhin [TB], which generalises Thouless’ 
formula to the case of correlated disorder. This discrepancy is a resonance 
effect and represents the “anomaly” in the presence of correlated disorder. 
On the other hand, disorder correlations modify the localisation length via 
the power-spectrum factor which is already present in the formula obtained 
by Izrailev and Krokhin in Ref. [TB]. By manipulating this factor, one can 
produce very strong localisation or effective delocalisation of the band-centre 
states. 

The modulation of the phase distribution close to the band centre is the 
hallmark of the anomaly that occurs there. Our analytical expressions show 
that positive, exponentially decaying correlations of the disorder flatten the 
invariant distribution, thereby reducing the anomaly, as observed in [T7]. The 
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opposite effect occurs when correlations decrease exponentially in magnitude 
but oscillate between positive and negative values. Such correlations strongly 
enhance the band-centre anomaly. To complete the picture, we consider a 
third type of correlations, which describe a lattice formed by two statistically 
independent and physically interpenetrating sublattices. We show that these 
correlations do not alter the modulation of the invariant distribution with 
respect to the case of uncorrelated disorder; however, they can increase or 
decrease the energy interval over which the anomaly is signihcant. 

This paper is organised as follows. In Sec. [2] we dehne the model under 
study and we introduce the Hamiltonian map approach. In Sec. |3] we show 
how to replace the random map for the phase variable with a correspond¬ 
ing stochastic differential equation. We proceed to derive general analytical 
expressions for the phase distribution and the localisation length in Sec. HI 
After recovering the known results for the limit case of uncorrelated disorder 
in Sec. O we consider the case of disorder with positive and exponentially 
decreasing correlations in Sec. |6l The case of correlations oscillating between 
positive and negative values and with exponentially decreasing magnitude is 
discussed in Sec. [71 A case of long-ranged correlations is analysed in Sec. [HI 
We hnally draw our conclusions in Sec. [H 

2 The Hamiltonian map approach 

2.1 Definition of the model 

We consider the ID Anderson model with weak and correlated disorder. The 
model is dehned by the Schrodinger equation 

V^n+l T 1 T (1) 

with random site energies We use energy units such that /2m = 1. We 
assume that 

{Sn) = 0 and {ei) = cx^. (2) 

We restrict our attention to the case of weak disorder, dehned by the condi¬ 
tion 

< 1. (3) 
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To complete the description of the statistical properties of the model ([I]), we 
introduce the normalised binary correlator 


{^n^n+l'} 


= x{i)- 


(4) 


Note that, in the weak-disorder limit, only the binary correlator (|1]) is re¬ 
quired to dehne the statistical properties of the disorder (unless one wants to 
go beyond the second-order approximation). We assume that the system is 
spatially homogeneous in the mean; for this reason the binary correlator (1T|) 
depends only on the distance I between the sites. We further assume disorder 
to be left-right symmetric on average so that x(/) is an even function of 1. 

We further assume that the binary correlator %(/) decreases quickly bey¬ 
ond a hnite length scale Ic- This condition was not used in previous second- 
order analyses of the Anderson model with weak disorder HD]. It is needed 
here, however, to avoid mathematical inconsistencies in the application of the 
special technique which is required to deal with the compound problem of 
the band-centre anomaly in the presence of correlated disorder (see Sec. H- 
We shall consider disorder to be weak enough that condition 


cr^/c “C 1 


(5) 


applies. In physical terms, this is equivalent to the assumption that the 
correlation length Ic is much shorter than the localisation length. To deal 
with the case of long-ranged correlations, we shall derive results valid for any 
hnite Ic and then consider the limit Ic —)■ oo, as discussed in Sec. [8l 


2.2 The Hamiltonian map 

The Hamiltonian map approach provides a useful way to study the structure 
of the electronic states of the Anderson model (JO]- The approach is based 
on the analogy between the quantum model ([T]) and a classical parametric 
oscillator with Hamiltonian 

^ ^ [1 + ^(t)] (6) 

where ^{t) is a succession of delta kicks of random strengths 

OO 

n=—oo 
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The integration of the dynamical equations of the oscillator (E]) over the 
period T between two kicks leads to the Hamiltonian map 

Xn+l = [cos (uT) - uiri Sm (cjT)] Xn + — siu (wT) Pn 

mu I ') 

Pn+l = —mu [sin (uT) + U^n cos (uT)] Xn + cos (uT) Pn- 

By eliminating the momenta from the map ([7]), one obtains the equation 

Xn+l + Xn-1 + U^n sin {uT) = 2 COS {uT) Xn, 

which coincides with the Schrodinger equation ([1]) for the Anderson model 
provided that 


E = 2 cos {uT) and Sn = u^n sin {uT). 

It is convenient to write the Hamiltonian map ([7]) in terms of the action-angle 
variables {Jn,0n), dehned by the equations 


Xn 

Pn 


2Jn 


sin On 


mu 


\/2muJn cos 9n- 


In this way one arrives at the map 

6n+l = 6n + uT + uin sin^ On + {u^nY siu^ On COS On 

+ o (a^) (mod 27r) 

Jn sin On cos On ~\~ u sin On'^ • 

In Eq. dH]) we have used the Landau symbol o(ct^) to denote neglected terms 
which, in the limit cr —?• 0, vanish faster than (see, e.g., [22])- In what 
follows we shall mostly omit the symbol o (cr^) and tacitly assume that the 
identities are correct within the limits of the second-order approximation in 
the disorder strength. 


2.3 The inverse localisation length 

The inverse localisation length (or Lyapunov exponent) is defined as 

A = lim — V In . 

N ^ '^n-l 
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Going to action-angle variables, one can write the Lyapunov exponent as 


A 


N 


lim — 

V—>oo N 


^In 


n=l 


+ lim ^ In 

Jyi—\ N—^oo iV 


sin 6 n 
sin 9 q 


(9) 


Except that at the band edge (where the angular variable tends to assume 
the values 0 and vr), the second term in the right-hand side (rhs) of Eq. ([9]) 
vanishes; one is therefore left with 

1 ^ J 

\= lim — In = (In (l - 2u^n sin On cos 9n + sin^ 6*^) )• (10) 

V^oo iV ^^ Jn-^ 
n=l ^ 

For weak disorder it is possible to expand the logarithm in the rhs of Eq. flTOj) 
and write 

2 

A = [1 - 2(cos(20„)) + (cos(40„))] - sin (20„)). (11) 

The noise-angle correlator (.^nsin {29n)) vanishes if the random site ener¬ 
gies are independent, but is different from zero for correlated disorder. It can 
be evaluated with the method presented in Ref. [19]; substituting the result 
in Eq. flTT]) one obtains 

2 

^ = ^-4r^{|l-2(cos(2»„)) + (cos(4«„))]H/(a,T) 

8sm (wT) 

+ [2(sin {29n)) - (sin (d^^))] Y (wT)} , 

where 

OD 

IE(a;) = 1-f 2 ^ x(/) cos (2a;/) (13) 

i=i 

is the power spectrum of the disorder and 

OD 

F(a;) = 2 ^ x(/) sin (2a;/) (14) 

i=i 

is the sine transform of the binary correlator fjTl) . 

To evaluate the averages of the trigonometric functions in the rhs of 
Eq. flT^ . it is necessary to determine the distribution of the angle variable 9. 
If uT does not lie too close to a value of the form vrp/g (with p and q integer 
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numbers), one can see from the map ([H]) that the angular variable quickly 
sweeps the [ 0 , 2tt] interval, thus producing a uniform invariant distribution, 

m - 

If this distribution is used to compute the averages of the trigonometric 
functions in Eq. flT^ . one immediately obtains the standard formula derived 
by Izrailev and Krokhin in Ref. |19] . 


Atk — 


a 


8 sin^ (wT) 


W{u:T). 


(15) 


When uT is a rational multiple of vr, however, the map (0) has almost periodic 
orbits whose existence manifests itself in the form of a periodic modulation 
of p{0). This is what happens at the band centre, which corresponds to the 
value uT = 7i/2. In order to obtain the correct localisation length close to 
the band centre, one must therefore determine the invariant distribution p{6) 
for the angular map ([ 8 ]) with uT ~ 7 r/ 2 . 


2.4 The angle map in a neighbourhood of the band 
centre 

We consider the case in which 

uT ~ with 5^0. (16) 

The corresponding energies lie close to the band centre, 

E = —2 sin 5 ~ —26. 

When the parameter uT takes the value flT 6 |) . the angular map ([ 8 ]) has almost- 
periodic orbits of period 4 and the difference 6*^+4 — is small. Iterating 
four times the angular map ([8]) leads to 

(^n +4 = 6n + 4:6 ---^ [1 — cos (20„)] H-i— - — [1 + cos {26n)\ 

2^2 ^ 

- Y [Ai + x(3)] sin {2dn) + ^ [2 - 3x(l) + 2^(2) - x(3)] sin (40„) 

(mod 27r). 


(17) 





In the rhs of Eq. (fT7|) we have neglected terms of order O (aS) and o (cr^). In 
order to determine the invariant measure p{0), it is useful to go to the con¬ 
tinuum limit and replace the random map ffT7|l with an appropriate stochastic 
differential equation. We devote the next section to the discussion of how 
this can be done. 


3 The continuum limit 

We need a systematic recipe to associate a stochastic differential equation to 
a random map of the form (IT7|) . Devising such a method is equivalent to End¬ 
ing an integration scheme for stochastic differential equations with coloured 
noise. Although there is an enormous literature on the numerical integration 
of stochastic differential equations with white noise (see |23] and references 
therein), much less is known on how to deal with differential equations with 
coloured noise. Equations with coloured noise are often reduced to systems 
of coupled equations with white noise with the standard trick of representing 
the coloured noise as the solution of an extra stochastic differential equation 
(see, e.g. [2ll |25]). In essence, this method replaces an equation of the form 

X = a{x) + b{x)7]{t), 

where ri{t) is a coloured noise, with the pair of coupled equations 

X = a(x) -f b{x)r](t) 

V = -7h(f) + Ht) 

with L{t) being a white noise. Numerical integration schemes have been 
created for equations of this kind (see, e.g., |2B]). This approach, however, 
can be applied only if the coloured noise is exponentially correlated, and we 
would like to avoid such a constraint. For this reason, we have derived a new 
integration scheme, which does not suffer from the same limitations. 

Let us consider a stochastic differential equation of the form 

N 

X = a{x,t) + '^^b^^\x,t)(i{t) (18) 

i=l 

where a{x,t) and the b^'^\x,t) are deterministic functions (with the functions 
f) being differentiable), while the (i{t) are stochastic processes with 
zero averages, 

(0(f)) = 0 (19) 
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and correlation matrix of the form 


{Ciit)Cjit + ^)) = Xij{r). (20) 

The matrix elements Xiji'i') Eq. fl20|) are assnmed to be even and decreasing 
fnnctions of the time difference r. Since we are interested in the case of 
weak noise, we will not specify the statistical properties of the (i{t) processes 
in fnrther detail. If the (i{t) are independent white noises, i.e., Ganssian 
“processes” with correlation fnnctions 

iCiiiXjii + r)) = 

we can interpret Eq. flTSl) as an informal way to write the Stratonovich 
stochastic differential eqnation 


N 

dx = a{x, t)dt + E t) o dWi- 

i=l 

Following [28] , we integrate the stochastic eqnation ffTSj) over a short time 
interval [t, t + e] and we obtain 




x{t + e) = x{t) + 


N 


a(x(r),r) + ^ 6W(x(r), r)Ci(r) 


2 = 1 


dr. 


A recnrsive application of this identity gives 




X 


(t + e) = x(t) + a(x(t),t)e + / dr Ci(T) 


db(d 




+ / dri I dr2 Ci(A)Ci(A) + • 




Neglecting the flnctnations of the noisy qnadratic term resnlts in the following 
integration scheme 


x(t + e) 


db^d 


x{t) + a{x{t),t)e + ^hj-^{x{t),t)b^^\x{t),t) 






/ Ci(rdr + ... 


( 21 ) 
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where the symbol Jjj stands for the integral 

hj = ^j (e - |r|) Xij(-|r|)dr. (22) 

It is now convenient to consider the discrete times tn = ne and to introdnce 
the compact notations 


Xn = X{tn), 


Qin ®^n)) 

db^n db^''‘'> , 

^nitn) • 


(23) 


dx dx 
We also define the new random variables 


= 


^n+1 


Ci(T)dT. 


(24) 


The statistical properties of the noises (i(t) define the corresponding prop¬ 
erties of the random variables fl2T)) . In particnlar, Eqs. flT^ and fl20l) imply 
that 

(Z®) = 0 

and 

= f (e - |r|) xp(r + fce)dr. 


We can now write Eq. fl?ID in the form of a map 


Xn-\-l Xj^ 


db[ 


(i) 




(25) 




The map fl2^ represents an integration scheme for the stochastic differential 
eqnation flTSl) . 

We now focns onr attention to the case in which the binary correlator (1201) 
has the form 

CXD 

{Ci{t)Cj{t + T)) = ^ Xij(4fc)5(r-fce), (26) 

k=—oo 

with Xij{4k) being a decreasing fnnction of the argnment 4k. In this case 
the correlation fnnction of the random variables fl24D becomes 

(Z® Z® J = eX.,(4k). (27) 
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Eq. 026 p also allows one to write the integrals 02 2 p in the form 

Substituting this identity in the map fl25ll and setting e = 1, we hnally obtain 
^n+i = Z«, (28) 

ij i 

The map fl28P contains the random variables Zn'^ with zero averages and cor¬ 
relation function dehned by Eq. fl27p . Our derivation shows that the map fESp 
represents an integration scheme for the stochastic differential equation flTKp . 
which contains noises Q{t) having vanishing averages and correlation func¬ 
tion fl26p . 

We can now address the inverse problem, i.e., how to associate a stochastic 
differential equation to a given a random map. Let us consider the map 

Xn+l = Xn + dn + (29) 


where a„ and bn'^ are short-hand notations dehned by Eq. fl2^ and the sym¬ 
bols Zn'^ represent random variables with zero average and binary correlator 
of the form (j2711 . The correspondence between the stochastic differential 
equation (ITsp and the map (j28ll implies that the random map (12^ can be 
read as an integration scheme for the stochastic differential equation 


1 _ 

X = a{x,t) - 

i 


(30) 


where the stochastic processes Q{t) have zero averages and correlation func¬ 
tion given by Eq. fl26l) . 

To test the correspondence between the map fl2^ and the stochastic dif¬ 
ferential equation fl30l) . we can apply it to the case of the random map (I8|) for 
the angle variable. In this case, the continuum limit is a differential equation 
of the form 

6 = uT + (^(t) siv? 6 (31) 
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where ({t) is a noise with zero average and correlation function 

CXD 

+ r)) {U^n+k)S (r - kT). (32) 

k=—oo 

Eqs. fl3Tl) and fl32|) dehne a dynamical system which was shown to be the 
continuous counterpart of the discrete Anderson model ([1]) for non-resonant 
values of the energy [2^. In Ref. |27], however, the correspondence between 
the two models was established only retrospectively, by computing separately 
the Lyapunov exponents of both systems and observing that they coincide. 
The integration scheme derived in this section, on the other hand, has en¬ 
abled us to predict a priori that the dynamical system is the continuous 
analogue of the discrete Anderson model ([1]). This is of crucial importance 
in the present case, where the continuum limit is instrumental in deriving an 
expression for the localisation length at the band centre. 


4 Invariant measure and localisation length 
in a neighbourhood of the band centre 

4.1 The stochastic differential equation 

Having established a correspondence between a generic random map of the 
form fl29l) and the stochastic differential equation fl30|l . we can now consider 
the continuum limit of the specihc map ffT7|l . As a hrst step, we can write 
the map flTTD in the form 

2 

On+A = 6'n -F 45 - y [x(l) -f- x(3)] sin {29^) 

- y [2 - 3x(l) + 2x(2) - x(3)] sin (40„) 

+ [1 - cos (29^)] + [1 + cos {29n)] 

with 

= gn +^n+2 ^^2) _ Sn+l+En+S ^ 


13 




The random variables fl34p have zero average and Eq. (jl]) implies that their 
correlation matrix has elements 

ZnUk) = = y [2x(4fc) + xiik + 2) + x(44 - 2)] 

(Z<‘' zZl^) = ^ [2x(4A^ + 1) + x(4A^ + 3) + X'(4«: - 1)] 

(Zf zZUk) = Y 12X(4*^ - 1) + X(4fc + 1) + X(4fc - 3)1. 

Following the approach discnssed in Sec. [3l we can interpret the map fl3^ 
as the integration scheme with time step 

e = 4T = 1 

of a stochastic differential eqnation of the form 

e = F^^\e) + F^^\e,t), ( 35 ) 

where {0) is the deterministic fnnction 

F(o)(0)=4<5-y [x(l) + x(3)]sin(20), 

while the stochastic part is given by 

F^^\6, t) = Ci(t) [1 - cos {26)] + C, 2 {t) [1 + cos {26)] . (36) 

In Eq. fl36|l . (i{t) and ( 2 {t) are two cross-correlated colonred noises with zero 
average and binary correlators 

OO 

{C(«)Cj(i + r)) = ^ {Z<-' zZ,t)S (t - k) (37) 

k=—oo 

with i,j = 1, 2. Note that the condition that the binary correlator (jTj) decays 
qnickly over distances larger than implies that the correlators 0371) also 
vanish over time scales r 3> Tc ~ Ic- 

We remark that, in the limit case of nncorrelated disorder, Eq. fl3^ re- 
dnces to 

d6 = 4:6dt + \ j~^ [1 ~ cos (26*)] o dWi -I- [1 + cos (26*)] o dfl^, 

which is the Stratonovich form of the ltd eqnation used in (21] to analyse the 
band-centre anomaly in the case of nncorrelated disorder. 
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4.2 Associated Fokker-Planck equation 

To proceed further, we make use of assumption ([5]), which we can now write 
in the equivalent form 

TeCT^ < 1. (38) 

As shown by Van Kampen [25], if condition fl38l) holds one can associate 
a Fokker-Planck equation to the stochastic differential equation ([35]). The 
statistical properties of the solution 9(t) of Eq. (I3^ can then be described 
in terms of a probability P{6,t) which obeys the associated Fokker-Planck 
equation 


|C,(»)P(«, ()| + ^ ICW)P(», ()l. (39) 

The drift and diffusion coefficients in Eq. f|3^ are given by 

poo ^Z7(l) 

Ci(e) = F<»)(«)+/ {—{e,t)F'-'He-\t-T))—dT (4o) 

v 0 

and 

poo 

C2ie)= / {F^^\e,t)F^^\e-\t-T))dT. (4i) 

Jo 

In Eqs. flTU]) and fITTD . the symbol 6^ represents the solution of the ordinary 
differential equation 

0t ^ 

with initial condition 6^ = 9. 

To simplify the mathematical expressions, in what follows we assume that 
Tc satishes not only condition fl38l) but also the additional condition 

< 1. (42) 

The present approach can be applied even if condition (1421) does not hold; 
assuming that it does, however, is mathematically convenient because the 
combination of conditions fl38p and (1421) ensures that 0* ~ 6* over time scales 
|f| < Tc and therefore one can replace 9~'^ with 9 in the integral expressions 
of the coefficients (140|) and (14T1) . This leads to a signihcant simplihcation of 
the mathematical expressions. From a physical point of view, condition (1421) 
is not very restrictive, because we are interested in the neighbourhood of the 
band centre, which corresponds to the limit 5 —)■ 0. 
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Carrying out the calculations, one obtains that the Fokker-Planck equa¬ 
tion fl39p takes the specihc form 


dP d 

- TW 


8>i-2W 


TT 


2hh(o) + iy 


TT 


sin(46') P{0,t) 

1 dP 

[1 + cos(4»)]j 1-, 


(43) 


where W is the power spectrum (IT^ and we have introduced the parameter 


X = 


26 




(44) 


4.3 The invariant measure 


We are interested in the stationary solution of the Fokker-Planck equa¬ 
tion (1431) . i.e., in the function 


p{9) = lim P{9, t), 


t^OO 


(46) 


which represents the invariant measure of the map ffT7|) . The distribution 
satishes the hrst-order differential equation 


dp 

d9 


(9) = A{9)p{9) + CB{9) 


(46) 


where C is an integration constant while 44(6*) and B{9) are the functions 

2W{'k/2) sin(40) - 8x 


and 


A{9) = 

B{9) = 


2fF(0) + hF(7r/2) [1 4- cos(40)] 
1 


2hF(0) + W(7r/2) [1 + cos(40)] ‘ 
The general solution of Eq. (146|) is 

^-8>cF{d) 


PiO) = 


^/2W{0) + hF(7r/2) [1 + cos(40)] 

c ‘ 


lo ^21T(0) + hF(7r/2) [1 + cos(40)] 


{\/2HV(0) + W'(ir/2)]p(0) 

d0 


(47) 
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The function F{6) in Eq. fITTj) is defined by the integral representation 


m = 


d(p. 


/o 21E(0) + lE(7r/2) [1 + cos(40)] 
Carrying out the integration one obtains the explicit expression 


F{9) = 


4^fE(0) [iy(0) + W{7r/2)] 


< arctan 


W{d) 


■ tan [29) 


W{0) + W{tt/2) 

(48) 

for (2n — l)7r/4 < 9 < (2n + l)7r/4 and n E Z. In Eq. fl48|) . arctan(a;) is a 
function with principal values in the interval [—7r/2,7r/2]. Note that 


+ Tin 


F{27i) = 


IT 


^W{0) [1E(0) + lE(7r/2)]' 


The constant C in Eq. fj47ll can be expressed in terms of p(0) with the 
help of the periodicity condition p{2tt) = p(0). The normalisation condition 
then determines the remaining integration constant p(0). In this way one 
arrives at the desired invariant measure 


p{9, h) = 


,-8xF(0) 


'•27r 




^21E(0) + lE(f) [1 + cos(40)] \}b ^2W{0) + lE(f) [1 + cos(40)] 


d(p 


+ e' 


87ri^(27r) 


^S>cF{4>) 


d0 


/o ywpyTW (| y [ rTSS ( 4 ^ j n(x) 


(49) 


where N{x) is a normalisation constant equal to 

^27r g-8xF(0) r ^27r 


38xF(</)) 


N{k) = 

_|_ g87rF(27r) 


0 ^j2W{d) + lT(f) [1 + cos(40)] \Je ^2W{d) + iy(f) [1 + cos(40)] 

fd pSxFCP) ] 

/ - drh I d9 

Jo y/2W{0) + lE(f) [1 + cos(40)] 

In Eq. (I49|) we have written p as a function of two arguments, 9 and x, to 
stress that the invariant distribution of the angular variable 9 depends on 
the energy parameter fHTll . 

A key feature of the invariant distribution is that it is F/2-periodic: 


dcj) 


P (^ + 1’^) = • 


(50) 


17 
























The demonstration of Eq. fISOp can be carried out along the same lines fol¬ 
lowed in the case of uncorrelated disorder [2T]. At first sight, Eq. (ISUl) may 
seem surprising because the deterministic part of the map fl33l) has stable 
zero-velocity points for 6 = 0 and 6 = ±7r. One could therefore expect prob¬ 
ability peaks to appear at these points and the distribution p{9, x) to be 
TT-periodic. However, a closer examination of the map (1551) reveals that the 
noisy terms are largest exactly for 6 = 0 and 6 = ±7r. The noise, therefore, 
scatters the angle variable away from the stable zero-velocity points. Since 
the deterministic and the noisy terms in the map (I55P carry the same weight 
at the band centre, the period of the invariant distribution ends up being 
7r/2 rather than tt. 

To conclude our analysis of the invariant distribution (I55p . we can consider 
its limit forms for x = 0 and x —)■ oo. At the exact band centre (x = 0), 
Eq. (Hop reduces to the signihcantly simpler form 


P{6, X = 0) = 


2K 


W{'k/2) 


1E(0) + lE(7r/2) 


U- 


2W{n/2) 


W{0) + lE(7r/2) 


[1 — cos(46')] 


(51) 

where the symbol K(fc) represents the complete elliptic integral of the hrst 
kind. In the opposite limit, i.e., for x oo, the asymptotic expansion of the 
general formula (H9l) gives 


p{9, x) 


1 


l + ^W 

4x 


(i) ^ ■ 


(52) 


Eq. (I52l) shows that, when the energy moves away from the band centre, the 
difference between the invariant distribution and its flat limit form falls off 
as the inverse of the energy. 


4.4 The localisation length 

Having determined the invariant measure fl49p . we can now turn our attention 
to the localisation length. Since p{6, x) is 7r/2-periodic, one has 

(cos(26*)) = (sin(26*)) = 0. 

Therefore, in a neighbourhood of the band centre, Eq. fll2p reduces to 


18 










with the functions W and Y dehned by Eqs. fll3p and fll4p . Taking into 
account that Y{'k/2) = 0 and neglecting terms of order one can 

write Eq. fl53P as 

A = y [1 + (cos (40n))] W , (54) 

where the average of the trigonometric function must be computed using the 
invariant measure flTOl) . 

Eq. flSrp represents the expression for the inverse localisation length in 
a neighbourhood of the band centre. It should be compared with the for¬ 
mula (IT^ for non-resonant energies which, close to the band centre, reduces 
to 

Aik = (I) . (55) 

As can be seen from Eqs. fl54p and 0551) . the modulation of the invariant 
distribution produces an anomaly in the inverse localisation length because 
the term (cos(46')) does not vanish. If the inverse localisation length fISTp is 
compared with the expression for uncorrelated disorder, it becomes obvious 
that spatial correlations have a twofold effect on the localisation length. On 
the one hand, the correlations manifest themselves via the power-spectrum 
factor iy(7r/2), which is present also in the standard formula (j55l) . On the 
other hand, the correlations modify the invariant distribution p(0, k) and 
therefore affect the value of the averaged cosine in Eq. fISTp . 

An explicit evaluation of the (cos(46')) is not possible in the general case, 
but it can be obtained at the exact band centre, where the invariant dis¬ 
tribution takes the simpler form fl5ip . In this case the inverse localisation 
length fl54p can be written as 

tr(^/2) ) 

M/(0) + iy(^/2) I 

- TVIO) ^ 

W-(>r/2) \ 

W'(O) + H2(?r/2) I 

(56) 

where K(fc) and E(/c) represent the complete elliptic integrals of the hrst and 
second kind. 

Eqs. (Hop and flSTl) constitute the central results of this paper. They 
provide general formulae for the invariant measure of the angular variable 
and for the localisation length of the Anderson model ([ip. To understand 


cr 


A(x = 0) = — < 


IE(0) + W 


TT 
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how the correlations of the disorder shape the localisation length and the 
invariant measnre, it is nsefnl to apply the general expressions (1491) and (1541) 
to specihc forms of correlated disorder. We devote the rest of this paper to 
this task. 


5 The case of uncorrelated disorder 


As a hrst application of the general formnlae obtained in Sec. HI we consider 
the limit case of nncorrelated disorder. When the site energies Sn are inde¬ 
pendent random variables, the binary correlator (H]) rednces to a Kronecker 
delta 

Xo (0 = ^10 

and the power spectrnm (IT^ takes a constant nnitary valne. The invariant 
distribntion p is obtained from the general expression (I49|) with W(0) = 
W(7r/2) = 1. One has 


p{e, x) = 


+ 


1 g-8xFo(0) 

N{k) .^3 + cos(46*) 



g8xFo(</>) 

'\/3T”cos(4^ 


d0 



g8>cFo (</>) 

-y/3 + cos(40) 


d0 


5 


(57) 


where N{x) is the normalisation factor and the fnnction Fq{6) is dehned by 
the integral expression 


j (68) 

Jq 3 + cos(40) 

Eq. (j57j) coincides with the resnlt originally obtained in |2T] . Taking into ac- 
connt that W(7r/2) = 1 for nncorrelated disorder, the general expression (jHTj) 
for the inverse localisation length rednces to 

2 

= y + (cos(40))], (59) 

where the average of the cosine fnnction mnst be compnted nsing the dis¬ 
tribntion (|571) . In Fig. [1] we compare the nnmerically compnted Lyapnnov 
exponent with the expression (15^ and with Thouless’ formula. 
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Figure 1: Inverse localisation length A versus in a neighbourhood of the 
band centre for uncorrelated disorder. The points correspond to the numer¬ 
ical data; the solid line to Eq. fl59|) : the dashed line to Thouless’ formula. 
The data were obtained for cx^ = 10“^. 


At the band centre, the invariant distribution (1511) and the inverse local¬ 
isation length (l56)) assume the well-known forms [20] 


p(e,^ = o) 


1 1 
K(l/\/2) -^/S -|- cos(46*) ’ 


(60) 


and 


A(x = 0) 





rr (3/4)1 

4 

K (l/x/2) 


Lr(i/4)J 


(61) 


In conclusion, the general expressions and flSTl) reproduce the correct 
results in the limit case of uncorrelated disorder. 
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6 Disorder with exponentially decaying, pos¬ 
itive correlations 


In this section and in those that follow we consider the Anderson model ([T]) 
with correlated site energies. Sequences of random variables {Sn} with arbit¬ 
rary binary correlations can be generated with the usual technique of hltering 
sequences of uncorrelated random variables |10j . 

We now focus our attention on the case of random site energies with pos¬ 
itive spatial correlations that decay exponentially with the distance between 
sites. The key result is that the band-centre anomaly is quickly suppressed 
for increasing values of the correlation length, as hrst observed in |T7]. As 
for the localisation length, it increases linearly with the correlation length, 
hoc oc Iq. 

We consider a binary correlator (0]) of the form 

Xi{l) = (62) 


The corresponding power spectrum (IT^ is 


OO 

l=—oo 


sinh(l//c) 

cosh(l//c) — cos(2a;) 


(63) 


The behaviour of the power spectrum (16^ is represented in Fig. [2] for various 
values of the correlation length. The values of the power spectrum (]63|1 at 
the boundaries of its domain are 




1 

W,{7r/2) 


sinh(l//c) 
cosh(l//c) — 1 


(64) 


Note that IFi(O) and IFi(7r/2) are, respectively, an increasing and a decreas¬ 
ing function of 

The invariant distribution for the angle variable is obtained by substitut¬ 
ing the values fl64p in the general expression 0491) . The analytical prediction 
matches well the numerical results, as can be seen in Figs. Eland 01 Fig. El 
shows how the resonance effect is progressively reduced as the energy moves 
away from the band centre. This effect is expected and is found also in the 
case of uncorrelated disorder [2T] . 

More surprising is the decrease of the modulation of the invariant meas¬ 
ure fl4^ with the correlation length. The suppression of the band-centre 


22 






Figure 2: Power spectrum fl6^ versus x for various values of the correlations 
length. 


anomaly as Ic is increased is particularly evident at the band-centre, where 
the invariant distribution (l5Til takes the form 

p(6*, x = 0) = — j — s - , (65) 

K [1-cos (40)] 

with 

For Ic —)■ 0, expression fl65l) differs from the corresponding formula (1601) for 
uncorrelated disorder by exponentially small terms of order As 

Ic increases, however, the modulation of p diminishes quickly, as shown by 
Fig. m To understand this effect, it is useful to observe that, for /c S> 1, 
expression fl65]) reduces to 


p{9,k = 0 ) 


1 



1 


cos (46^) -I-... 


( 66 ) 
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Figure 3: Invariant distribution p{d, k) for various values of k. The disorder 
exhibits positive exponential correlations of the form fl62ll with lc = 2. Lines 
correspond to the theoretical formula fH9l) : points to the numerical data. The 
numerical data were obtained for cx^ = 10“^. 


The asymptotic form (l66|) matches relatively well the numerical data already 
for values of /c ^ 10, as shown by Fig. [5l 

The inverse localisation length can be computed with the help of Eqs. (H9H 
and (I5TD . supplemented by the specihc values (IMD of the power spectrum. In 
Fig. Owe compare the numerically computed Lyapunov exponent with the 
theoretical expression fl54p and with the formula fllSp obtained by Izrailev 
and Krokhin. Note that the differences between the numerical values of A 
and the values predicted by Eq. flSTp are of order ~ 10“®, well within the 
error intrinsic to the second-order approximation used in the theoret¬ 
ical calculations. 

At the exact band centre, a relatively simple analytical expression for the 
inverse localisation length can be obtained by substituting the values (IMP in 
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Figure 4: Invariant distribution p{d, x = 0) at the band centre for various 
values of the correlation length Ic- The disorder exhibits positive exponen¬ 
tial correlations of the form fl6^ . The lines correspond to the theoretical 
formula ([65]); the points to the numerical data. The data were obtained for 

= 10-T 


Eq. (I56|) . One obtains 
A(x = 0) = 


a 


4sinh (l//c) 


2 cosh (1/Zc 


e ( V »’-(«/2 




— cosh (l//c) — 1 


Note that, for /c ^ 1, the inverse localisation length (1671) becomes 


(67) 


A(x = 0) ~ 


a 


Therefore the extension of the band-centre state increases linearly with the 
correlation length Ic- 


25 











0.1596 


P 


0.1595 

0.1594 

0.1593 

0.1592 

0.1591 

0.159 

0.1589 

0.1588 


Numerical data 
Analytical formula 



-71 -37r/4 -7r/2 -7r/4 0 7r/4 7r/2 37r/4 k 

e 


Figure 5: Invariant distribution p(P, x = 0) at the band centre for disorder 
with correlations of the form fl62ll with Ic = 10. The line corresponds to the 
asymptotic formula fl66l) : the points to the numerical data. The data were 
obtained for = 10“^. 


7 Exponentially decaying correlations with os¬ 
cillating sign 


We now consider correlations which decay exponentially with the distance 
between sites, but whose sign oscillates. Contrary to the previous case, the 
anomaly is now reinforced as the correlation length increases. 

Mathematically, the correlations of the site energies have the form 

x,(l) = (68) 


The corresponding power spectrum is 


OO 

in(x) = 

l=—oo 


sinh(l//c) 

cosh(l//c) + cos(2a:) 


(69) 
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Figure 6: Inverse localisation length A versus in a neighbourhood of the 
band centre for disorder with correlations of the form fl62l) and /c = 1- The 
points correspond to the numerical data; the solid line to Eq. fl5T|) : the dashed 
line to the formula ([15]). The numerical data were obtained for = 10“^. 


The behaviour of the power spectrum fl69|) is represented in Fig. [T] Eq. fl6^ 
implies that 


hF2(0) 


1 

fF2(7r/2) 


sinh(l//c) 
cosh(l//c) + 1 


(70) 


Note that 

W 2 (x) = VFi (I - x^ 


and that, therefore, the formulae for the invariant distribution and the inverse 
localisation length for this case can be obtained from the expressions for the 
previous case with the exchanges IFi(O) O W2{t^/2) and Wi{t:/2) o 114(0). 

As before, the invariant distribution for the angle variable is obtained 
by substituting the specific values fl7U]) of 1F(0) and lF(7r/2) in the general 
expression fl49p . The resulting analytical formula is corroborated by the 
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Figure 7: Power spectrum fl6^ versus x for various values of the correlations 
length. 


numerical results, as shown by Figs. [8]and[9l Fig. |8] shows how correlation 
of the form fl68|) produce a strong modulation of the invariant measure at the 
band centre which is gradually reduced as the energy moves away from the 
band centre. 

Fig. m on the other hand, shows how the invariant distribution develops 
conspicuous peaks as Ic is increased. This behaviour is entirely consistent 
with the theoretical predictions. At the band centre, in fact, the general 
expression fl5ip for the invariant distribution takes the form 

p(0,x = O) = — . ^ ^ (71) 

K (^Vv?+(/e)/2) V4-(p+(y [1-cos (40)] 

with 

^ ^ + COSh(l/iJ' 
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Figure 8: Invariant distribution p{9, x) for various values of the parameter 
X. The lines correspond to the theoretical formula (09]); the points to the 
numerical data. The data were obtained for disorder with correlation of the 
form fl68|) with = 5. The disorder strength in numerical calculations was 
CT^ = 10-3. 


The asymptotic behaviour of the distribution fITT]) for ^ 1 is 

p{6, X = 0) ~ — 






(72) 


Note that for 


^ TT TT 

6^ ~ —h n— with n G Z 

4 2 


the distribution fl72|) assumes the value 

/TT 


















Figure 9: Invariant distribution p{d, x = 0) at the band centre for various 
values of the correlation length Ic. The disorder exhibits correlations of the 
form fl68|) . The lines correspond to the theoretical formula fl7T]) : the points 
to the numerical data. The data were obtained for = 10“^. 


which diverges for 3> 1. This entails that the distribution flTT]) develops 
four sharp maxima as Ic is increased, in agreement with the data of Fig. |9l 

After inserting the values fl70|) in the expression for the invariant 
distribution, one can evaluate the rhs of Eq. fl54p and obtain the inverse loc¬ 
alisation length. The results agree with the numerical data, as shown by 
Fig. [ini The slight discrepancy between numerical data and the theoretical 
expression is of order O(cr^) and should probably be attributed to the neg¬ 
lected fourth-order correction in Eq. fl5T|) . Fig. [10] graphically shows that 
correlations of the form (l68ll increase the relative deviation of the Lyapunov 
exponent from the value predicted by the formula (lT5ll obtained by Izrailev 
and Krokhin. This is consistent with the very large peaks that the invariant 
distribution develops in the present case. 

Correlations of the form fl68|) have also the effect of enhancing the local- 
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Figure 10; Inverse localisation length A versus in a neighbourhood of the 
band centre for disorder with correlations of the form fl68ll with Ic = 5. The 
points correspond to the numerical data; the solid line to Eq. fl5T|) : the dashed 
line to the formula ([15]). The numerical data were obtained for = 10“^. 


isation of the electronic states for increasing values of the correlation length. 
This is a consequence of the fact that the factor lF(7r/2), dehned by Eq. (1701) . 
is an increasing function of Ic- This conclusion is conhrmed by the explicit 
formula for the inverse localisation length at the exact band centre. For 
E = 0, Eq. fl56l) becomes 


A(x = 0) = 


a 


4sinh {1/lc 


2 cosh (l//c) 


E (a/ (p+{lc)/‘2 


k(^cp4Q/2 


— cosh (l//c) + 1 


(73) 

As in the previous case, in the limit Ic ^ 0 Eq. fl73l) differs from its counter¬ 
part fl6Tl) only by vanishing terms of order In the limit Ic 3> 1, on 
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the other hand, the inverse localisation length fl73p diverges as 


A(x = 0) 


2/c 1 

T [ln(8/c) ~ ^ ^ ■ 


In this case the band-centre state becomes strongly localised as increases. 


8 The case of a composite lattice 


We now consider the case of a ID chain with random energies whose binary 
correlator satishes the condition 


X(2/ + 1) = 0 for / e Z. (74) 

A seqnence {Sn} with correlations of the form fl74)) can be obtained by mixing 
two independent random seqnences {«„} and {/?„} with the same statistical 
properties. More precisely, one assnmes that 

(an) = {/3n) = 0 

and that the binary correlators are 

(anan+i) = Wnf3n+i) = ^nd {an(3m) = 0, 

where x(/) is an arbitrary fnnction, satisfying the condition that is qnickly 
decays for I ^ Ic- One can then dehne the random energies by setting 

S2n and £2n-\-l 


In physical terms, the chain is split in two independent and interpenetrating 
snblattices. 

Condition fl74|) implies that the power spectrum must satisfy the identity 


H-(0) = W . 


As a consequence, the invariant distribution (1491) assumes the particularly 
simple form 


p{9, x) 


1 g-8xFo(0) 

N{x) ^3 + cos(46') 



p8xFo(<j>) 

, d0 

^/3 + cos (40) 



p8kFoi(t>) 

, d0 

-^3 + cos(40) 




(75) 


32 













where iV is a normalisation constant, Fq{9) is defined by Eq. fISS]) . and 


X 

1E(0)' 


(76) 


The distribution fl75|) has the same form of the distribution fl57|) for uncorrel¬ 
ated disorder, the only difference being that the parameter of the latter 
is replaced by the rescaled parameter fl76p in the former. In physical terms, 
this means that the correlations do not modify the modulation of p, but they 
alter the scale over which the distance of the energy from the band centre 
is measured. This entails that the anomaly extends over a larger energy 
interval if 1T(0) > 1 and is restricted to a shrunken region if 1T(0) < 1. 

At the exact band centre, x = x = 0 and the invariant distribution 
reduces to the form flBUjl . Correspondingly, the inverse localisation length 
becomes 


X{x = 0 ) 


a 


4 


„ E(1/72) 

K (1/V2) 



(77) 


8.1 Long-ranged correlations 

We now focus our attention on a specific type of correlator fulfilling condi¬ 
tion fl74|l . i.e., 

1 -I- (—1)^ sin (2aZ) 

2 2a/ 

where the parameter a lies in the interval [0,7r/4]. The correlator fl78|) does 
not decreases quickly for I > 1^, and therefore does not satisfy one of the 
conditions used to derive our analytical results. It can be fitted in our the¬ 
oretical framework, however, if it considered as the limit form for /c —)■ cxo of 
the correlator 

_ 1 + (-1)^ sin (2a/) 

2 2a/ ^ Wc y ■ 

The power spectrum of this correlator is 






W^{x) 


+ 



arctan 


< arctan 

sin (2a -|- 2x) 


- arctan 

sin (2a -|- 2x) 

gl/ic _ cos (X 

la -|- 2x) 


gl/Zc _|_ cos (z 

la -|- 2x) 

sin (2a — 2x) 


sin (2a — 2x) 

\ 

qI/Ic — cos (2a — 2x) 


giAc cos (2a — 2x) 

1 





(80) 
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In the limit Ic ^ oo the power spectrum (ISUD tends to the form 


Wsix) 


71 

4a 

0 


if a: G [0, a] U 

TT 


if X E 


71 71' 

- a, — 

2 2 . 


’ 2 


a 


(81) 


The behaviour of the power spectrum flHOl) with a = tt/IO is represented in 
Fig. [TT] for various values of Ic- In the limit Ic —>■ oo, the power spectrum 



V 


Figure 11: Power spectrum fl80|) with a = tt/IO for various values of the 
correlations length. For /c —t oo, mobility edges appear at Xi = tt/IO and 
X2 = 2tt/5. 


vanishes for Xi = tt/IO < x < X 2 = 27t/2; according to the standard for¬ 
mula ([T5D, this generates mobility edges at Ei ~ ±0.618 and E 2 — ±1.902. 

When the binary correlator takes the form 0791) . the invariant distribu¬ 
tion is given by the expression fl75|l . This is true for any hnite value of the 
correlation length; under the reasonable assumption that p{9, x) should be 
a continuous function of Ic, one can conclude that the invariant distribution 
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keeps the form (175]) even in the limit —)■ oo, i.e., when the binary correl¬ 
ator is given by Eq. fl78p . The numerical data corroborate this conclusion, 
as shown by Fig. [121 As can be seen, the numerical data match well the 
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Figure 12: Invariant distribution p{9, k) for disorder with correlations of the 
form (1751) with a = vr/lO. The data correspond to the values x = 0 (band 
centre) and k = A. Lines correspond to the theoretical formula (I75|) : points to 
the numerical data. Also represented are the numerical data for the invariant 
distributions for uncorrelated disorder and values x = 0 and x = 4 of the 
parameter (i44|) . In numerical computations the disorder strength was set at 
ct2 = 10-2. 

theoretical distribution (I75|) both at the band centre (x = 0) and away from 
it (x = 4). In Fig. [121 we also plot the numerically obtained invariant distri¬ 
butions for uncorrelated disorder with x = 0 (band centre) and x = 4. The 
data show that, as expected, these two distributions collapse on the distribu¬ 
tions for correlated disorder with x = 0 and x = 4. Note that in the present 
case value IF(0) = 5/2 so that x = 10 when x = 4. 

We can now consider the behaviour of the localisation length. In Fig. [131 


X - X - Q, Corr. Dis. (Analytical) - 

X - X - 0, Corr. Dis. • 

X -0, Uncorr. Dis. t 

X - A {x = 10), Corr. Dis. (Analytical) . 

X - A{x = 10), Corr. Dis. * 

X - A, Uncorr. Dis. * 
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we plot the Lyapunov exponent as a function of the energy. One can eas¬ 
ily see that the formula 0151) proposed by Izrailev and Krokhin agrees well 
with the numerical data; in particular, effective mobility edges arise where 
expected. Discrepancies between the numerical data and the expression flT5|l 
appear at the band centre and at the band edges, however, where anomalies 
occur. The band-centre anomaly is represented in greater detail in Fig. [T31 



-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 

E 


Figure 13: Inverse localisation length A versus E for disorder with correlations 
of the form 0781) . The points correspond to numerical data; the dashed line 
to the standard formula ([15]). The data were obtained for cr^ = 10“^. 


where we compare the numerical data with the standard formula OT^ and 
with Eq. 05^ . The use of the expression 05^ . rather then OMD . is due to the 
fact that in this case the anomaly extends, as predicted, over a larger energy 
interval and, therefore, the term cos^((5) in Eq. 0531) cannot be approximated 
with unity as done so far. Fig. ITTl conhrms that our theoretical results work 
rather well even for long-ranged correlations of the form OTHl) . The exten¬ 
sion of the energy interval with anomalous behaviour can be clearly seen if 
one compares the data represented in Fig. HH with those corresponding to 
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Figure 14: Inverse localisation length A versus in a neighbourhood of 
the band centre for disorder with correlations of the form fl78l) . The points 
correspond to the numerical data; the solid line to Eq. fl53|) : the dashed line 
to the standard formula ([15]). The data were obtained for cx^ = 10“^. 


uncorrelated disorder shown in Fig. [H 

To conclude the discussion of long-ranged correlations, we can consider 
the case of disorder with correlations of the form 


xS) 


sin (2a/) 
(2a/) 




The limit form for /c —)■ cxo is 


X4(/) 


sin (2a/) 
(2a/) 


cos 



(82) 
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The power spectrum corresponding to the correlator 0821) is 


WJx) = 1 + — < arctan 
4a [ 

arctan 


gi/^c ^2x — 2a) 

cos {2x — 2a) 

+ sin {2x + 2a) 


+ arctan 




sin {2x + 2a) 


cos {2x + 2a) 
In the limit Ic ^ oo the power spectrum 


arctan 


^/Ic _ 


cos (2x + 2a) 
sin {2x — 2a) 


cos (2x — 2a) 


(83) 


tends to the form 


TT 


Wi{x) = 


A if 2: G 
TT — 4a 

0 if X e 
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-— a, — + a 
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“ TT 
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TT TT 
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We thus obtain a power spectrum which is complementary with respect to 
the case described by Eq. flSTl) . The behaviour of the power spectrum fl8^ 
with a = vr/lO is represented in Fig. [15] for various values of Ic- 

As can be seen from Eq. flH^ . for increasing values of Ic the power spec¬ 
trum tends to zero at the boundaries of the domain [0,7r/2]. This entails 
that the anomaly at the band centre becomes unstable as Ic ^ oo. In 
fact, although at the exact band centre the invariant distribution keeps the 
form (jbOj) . the fact that hE(0) = hE(7r/2) —)■ 0 implies that the rescaled para¬ 
meter (i76|) grows very quickly for any infinitesimal deviation of the energy 
from the band centre. Therefore the invariant distribution becomes uniform 
very fast when the energy moves away from the band centre. 

We can conclude that, depending on the value of the power spectrum at 
the band centre, correlations satisfying condition fl74p can either strengthen 
or weaken the band centre anomaly. They do not enhance or suppress the 
modulation of the invariant distribution, but they can widen or shrink the 
neighbourhood of the band centre where the invariant distribution is signi¬ 
ficantly non-uniform. 


9 Conclusions 

In this paper we have studied the band-centre anomaly in the ID Anderson 
model with weak correlated disorder. Our perturbative analysis used two es¬ 
sential tools: the Hamiltonian map approach and the continuum limit. The 
Hamiltonian map approach interprets the spatial structure of the electronic 
states in terms of the time evolution of a classical parametric oscillator. The 
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Figure 15: Power spectrum fl8^ with a = vr/lO for various values of the 
correlations length. As Ic —)■ oo, mobility edges arise at xi = tt/IO and 
X2 = 27r/5. 


dynamical evolution of the angle variable of this oscillator is dictated by the 
random map flTTD . Replacing this map with a corresponding stochastic differ¬ 
ential equation is a crucial step that allowed us to derive our analytical results 
and that required the elaboration of the specihc integration scheme fl28l) . We 
obtained analytical expressions for the invariant distribution of the phase 
variable and for the localisation length. These results are valid for weak dis¬ 
order with arbitrary correlations and generalise the formulae obtained in [21] 
for the case of uncorrelated disorder. 

When disorder is uncorrelated, the invariant distribution of the phase 
variable, which is uniform for non-resonant values of the energy, becomes 
modulated for energies lying in a neighbourhood of the band centre. This 
modulation, in turn, generates a deviation of the inverse localisation length 
from the values predicted by Thouless’ formula. In qualitative terms, this 
picture holds also when disorder displays spatial correlations. From a quant- 
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itative point of view, however, the size and extension of the resonance effect 
can be dramatically altered. Two extreme cases are discussed in Sec. O and [71 
In the hrst case, the site energies exhibit positive correlations which decay 
exponentially with the distance between sites. In this case the resonance 
effect is suppressed upon increasing the correlation length Ic- In the second 
case, the correlations between site energies also decrease exponentially in 
magnitude, but oscillate between positive and negative values. In this case, 
increasing the correlation length Ic strongly enhances the modulation of the 
invariant measure, which tends to a sum of four delta peaks for /c 3> 1. 
Correspondingly, the difference between the inverse localisation length and 
the value predicted by the formula derived by Izrailev and Krokhin increases 
with Ic- The specific long-range correlations analysed in Sec. jS] do not alter 
the modulation of the invariant distribution with respect to the case of un¬ 
correlated disorder, but can strengthen or weaken the anomaly in a different 
way, i.e., they can enlarge or restrict the interval of the energy in which the 
resonant effect is relevant. 

In conclusion, we have shown that correlations of the disorder can alter 
the band-centre anomaly very strongly and in a variety of ways. In particular, 
specihc correlations can suppress or magnify the resonance effect at the centre 
of the energy band. 
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